Metadynamic sampling of the free energy landscapes of proteins coupled with a Monte Carlo algorithm

Metadynamics is a powerful computational tool to obtain the free energy landscape of complex systems. The Monte Carlo algorithm has proven useful to calculate thermodynamic quantities associated with simplified models of proteins, and thus to gain an ever-increasing understanding on the general principles underlying the mechanism of protein folding. We show that it is possible to couple metadynamics and Monte Carlo algorithms to obtain the free energy of model proteins in a way which is computationally very economical.Comment: Submitted to Gen

Room Temperature Organic Superconductor?

The electron--phonon coupling in fullerene C28 has been calculated from first principles. The value of the associated coupling constant lambda/N(0) is found to be a factor three larger than that associated with C60. Assuming similar values of the density of levels at the Fermi surface N(0) and of the Coulomb pseudopotential for C28-based solids as those associated with alkali-doped fullerides A3C60, one obtains Tc(C28) \approx 8 Tc(C60).Comment: 10 pages, 2 figure

A folding inhibitor of the HIV-1 Protease

Being the HIV-1 Protease (HIV-1-PR) an essential enzyme in the viral life cycle, its inhibition can control AIDS. The folding of single domain proteins, like each of the monomers forming the HIV-1-PR homodimer, is controlled by local elementary structures (LES, folding units stabilized by strongly interacting, highly conserved, as a rule hydrophobic, amino acids). These LES have evolved over myriad of generations to recognize and strongly attract each other, so as to make the protein fold fast and be stable in its native conformation. Consequently, peptides displaying a sequence identical to those segments of the monomers associated with LES are expected to act as competitive inhibitors and thus destabilize the native structure of the enzyme. These inhibitors are unlikely to lead to escape mutants as they bind to the protease monomers through highly conserved amino acids which play an essential role in the folding process. The properties of one of the most promising inhibitors of the folding of the HIV-1-PR monomers found among these peptides is demonstrated with the help of spectrophotometric assays and CD spectroscopy

Molecular details of dimerization kinetics reveal negligible populations of transient µ-opioid receptor homodimers at physiological concentrations.

Various experimental and computational techniques have been employed over the past decade to provide structural and thermodynamic insights into G Protein-Coupled Receptor (GPCR) dimerization. Here, we use multiple microsecond-long, coarse-grained, biased and unbiased molecular dynamics simulations (a total of ~4 milliseconds) combined with multi-ensemble Markov state models to elucidate the kinetics of homodimerization of a prototypic GPCR, the µ-opioid receptor (MOR), embedded in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/cholesterol lipid bilayer. Analysis of these computations identifies kinetically distinct macrostates comprising several different short-lived dimeric configurations of either inactive or activated MOR. Calculated kinetic rates and fractions of dimers at different MOR concentrations suggest a negligible population of MOR homodimers at physiological concentrations, which is supported by acceptor photobleaching fluorescence resonance energy transfer (FRET) experiments. This study provides a rigorous, quantitative explanation for some conflicting experimental data on GPCR oligomerization

Sirolimus Therapy for Patients With Lymphangioleiomyomatosis Leads to Loss of Chylous Ascites and Circulating LAM Cells

A young woman received a diagnosis of abdominal, sporadic lymphangioleiomyomatosis (LAM) and multiple abdominal lymphangioleiomyomas and was referred for recurrent chylous ascites responding only to a fat-free diet. On admission, pulmonary function test (PFT) results showed a moderate reduction in the transfer factor for carbon monoxide with normal exercise performance. The serum vascular endothelial growth factor D (VEGF-D) level was 2,209 pg/mL. DNA sequences, amplified at loci kg8, D16S3395, D16S3024, D16S521, and D16S291 on chromosome 16p13.3, showed a loss of heterozygosity (LOH) only for kg8. Fat-free total parenteral nutrition in association with sirolimus (2 mg po daily) was initiated. Serum sirolimus levels were maintained at concentrations between 5 and 15 ng/mL. After 1 month, reintroduction of a low-fat oral feeding was achieved without recurrence of ascites. PFT results were stable. Interestingly, clinical improvement was associated with a reduction in the VEGF-D serum level (1,558 pg/mL). LOH at the kg8 biomarker in blood LAM cells was no longer detected

Deriving amino acid contact potentials from their frequencies of occurence in proteins: a lattice model study

The possibility of deriving the contact potentials between amino acids from their frequencies of occurence in proteins is discussed in evolutionary terms. This approach allows the use of traditional thermodynamics to describe such frequencies and, consequently, to develop a strategy to include in the calculations correlations due to the spatial proximity of the amino acids and to their overall tendency of being conserved in proteins. Making use of a lattice model to describe protein chains and defining a "true" potential, we test these strategies by selecting a database of folding model sequences, deriving the contact potentials from such sequences and comparing them with the "true" potential. Taking into account correlations allows for a markedly better prediction of the interaction potentials

Designability of lattice model heteropolymers

Protein folds are highly designable, in the sense that many sequences fold to the same conformation. In the present work we derive an expression for the designability in a 20 letter lattice model of proteins which, relying only on the Central Limit Theorem, has a generality which goes beyond the simple model used in its derivation. This expression displays an exponential dependence on the energy of the optimal sequence folding on the given conformation measured with respect to the lowest energy of the conformational dissimilar structures, energy difference which constitutes the only parameter controlling designability. Accordingly, the designability of a native conformation is intimately connected to the stability of the sequences folding to them.Comment: in press on Phys. Rev.

Simple models of protein folding and of non--conventional drug design

While all the information required for the folding of a protein is contained in its amino acid sequence, one has not yet learned how to extract this information to predict the three--dimensional, biologically active, native conformation of a protein whose sequence is known. Using insight obtained from simple model simulations of the folding of proteins, in particular of the fact that this phenomenon is essentially controlled by conserved (native) contacts among (few) strongly interacting ("hot"), as a rule hydrophobic, amino acids, which also stabilize local elementary structures (LES, hidden, incipient secondary structures like $\alpha$--helices and $\beta$--sheets) formed early in the folding process and leading to the postcritical folding nucleus (i.e., the minimum set of native contacts which bring the system pass beyond the highest free--energy barrier found in the whole folding process) it is possible to work out a succesful strategy for reading the native structure of designed proteins from the knowledge of only their amino acid sequence and of the contact energies among the amino acids. Because LES have undergone millions of years of evolution to selectively dock to their complementary structures, small peptides made out of the same amino acids as the LES are expected to selectively attach to the newly expressed (unfolded) protein and inhibit its folding, or to the native (fluctuating) native conformation and denaturate it. These peptides, or their mimetic molecules, can thus be used as effective non--conventional drugs to those already existing (and directed at neutralizing the active site of enzymes), displaying the advantage of not suffering from the uprise of resistance

Electron-phonon interaction in C70

The matrix elements of the deformation potential of C$_{70}$ are calculated by means of a simple, yet accurate solution of the electron-phonon coupling problem in fullerenes, based on a parametrization of the ground state electronic density of the system in terms of $sp^{2+x}$ hybridized orbitals. The value of the calculated dimensionless total electron-phonon coupling constant is $\lambda\approx0.1$, an order of magnitude smaller than in C$_{60}$, consistent with the lack of a superconducting phase transition in C$_{70}$A$_3$ fullerite, and in overall agreement with measurements of the broadening of Raman peaks in C$_{70}$K$_4$. We also calculate the photoemission cross section of C$_{70}^-$, which is found to display less structure than that associated with C$_{60}^-$, in overall agreement with the experimental findings.Comment: To be published in Phys. Rev.